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Introduction 


Modeling of chemistry in Computational Fluid Dynamics can be 
-he most time-consuming aspect of many applications. If the entire 
set of elementary reactions is to be solved, a set of stiff 
ordinary differential equations must be integrated. Some of the 
reactions take place at very high rates, requiring short time 
steps, while others take place more slowly and make little progress 
in the short time step integration. 

Historically, the problem has been approached in several ways: 

1) Single Step - Instantaneous Reaction: While computationally 
simple, this technique will over-predict conversion since the 
equilibrium point will be surpassed. 

2) Total Equilibrium: This is the assumption used in ODE and TDE 
at MSFC. The assumption may be adequate in describing the overall 
performance of an engine, but may fail to provide the detail 
required for the spatial resolution which CFD analysis is to 
provide. That is, the assumption may not be valid for the local 
conditions resulting from the CFD calculations. 

3) Reduced Mechanism - Finite Rate: The choice of the appropriate 
reduced set is difficult, since the local conditions are not 
constant throughout the CFD calculations. Since the importance of 
the elementary reactions may change within the computational space, 
the global set chosen may not be appropriate for the local 
conditions. 

A) Partial Equilibrium - Finite Rate: A portion of the reactants 
are assumed to be in equilibrium, while the remainder are 
integrated in finite rate kinetics. This approach removes the very 
short time step calculations and allows integration of fewer 
equations at longer time steps. Since there is an interaction 
between the equilibrated and non-equilibrated species, an 
equilibrium calculation is performed at each time step. KIVA (Los 
Alamos) utilizes this approach. However, the proper choice of 
equilibrated and non-equilibrated reactions may change with the 
local conditions. 

The goal of this work is to develop a procedure to 
automatically obtain sets of finite rate equations, consistent with 
a partial equilibrium assumption, from an elementary set 
appropriate to local conditions. The sets can be applied to the 
appropriate regions within the CFD space where the total 
equilibrium assumption is inappropriate. 
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Figure 1. Computational grid used for the FI engine 
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Since the variation in mixture ratio needed to be modeled for 
testing of the reaction reduction scheme, an approximation of the 
complicated injector element location (Figure 2.) was included in 
the input file. Since the LOx had a 20 degree impingement angle 
and the RP-1 had a 15 degree impingement angle, the droplets were 
introduced into cells 2 and 3 respectively. Also the variation in 
or ifice are per ring was also modeled to a certain extent, although 
the grid was not dense enough for a one-to-one correspondence with 
the injector rings. The presence of the baffles was ignored. A 3D 
solution would be required to take the presence of the radial 
baffles into account. 


Injector Pettart*: 

dndXat 

. full Cl wane type ^ 

• i^i ipMu a*t 

15 

7.H 

- Orifice Bloater, m 

10. ♦ 

- O-lftc. Bpacln*. m 


. fetdtier Cl*-** a 

20 

. i_i i «*lf-An*le. 

4.15 

- trlfice OtMotar, m 

- cpRf Ice lpectr«. m 

10. « 


17.1 

. Injection artteltr. 

455 

• fuel-tide Prmmn *rap, j 

. teldlter injection Veteclty. 

40.) 

n»i 

• fraction ^ t 

S.2 

70 

. f recti** ef fetal #,r * n **‘ 1 


ttebUlty *■*■*•* 

11 

• a \M»r *f Teete 

* Average Be^ T, »* * 

a 

?*oi-actari.Mc Velocity efficiency. 1 

no 



* o!jac«»t to radial befflee -ora ».2 

caa directly 

«rd *.» — 




Axial faa d bote* to *9, *15. *23. *27. *31. 
oxkfizar rinQ* r**trictad. 


Figure 2. Injector configuration used in the CFD solution. (From 
Oefelein and Yang, 1992) 


The inclusion of the variation of mixture ratio resulted in 
variation within the solution necessary to test the reaction 
reduction scheme. (Figure 3) However the actual test of the 
reduction procedure has yet to be performed due to time 
constraints . 
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Figure 3. Temperature predictions near the injector. 

4 • Conclusions and Future Work 

A two-phase solution of the FL engine with equilibrium in the 
vapor phase was found which showed enough variation wiThin the 
flald , to expect reduced reaction sets priced with the 

computation'll 0 grid . out * lned («***». 1991, to" vary with the 

of the he Mn°aT cp n D Cl “^°" - has ? et l? be confirmed by implementation 
_ , . . lna i D solution in the reaction reduction scheme 

selection of representative cells with the predicted inllt 
conditions is the next step in furthering this work. 

A two-phase CFD code (Refleqs, CFDRC) can provide new insiaht 
as to the combustion within the F-l engine. Better modeling of the 
impinging injectors should yield more" meaningful results. 
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